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ABSTRACT 

Flux-dominated solar dynamo models, which have demonstrated to be quite successful in reproducing most of the observed features of the large 
scale solar magnetic cycle, generally produce an inappropriate latitudinal distribution of the toroidal magnetic fields, showing fields of large 
magnitude in polar regions where the radial shear has a maximum amplitude. Employing a kinematic solar dynamo model, we here explore the 
contribution of both the radial and the latitudinal shear in the generation of the toroidal magnetic fields by varying the shape and the thickness 
of the solar tachocline. We also explore the effects of the diffusivity profile of the convective zone. Considering the shear term of the dynamo 
equation, (B p ■ V)Q=B r dQ./dr + Bo/rdQ./d9, we find that the latitudinal component is always dominant over the radial component at producing 
toroidal field amplification. These results are very sensitive to the adopted diffusivity profile, specially in the inner convection zone (which is 
caracterized by the diffusivity r] c and the radius r c of transition between a weak and a strong turbulent region). A diagram of the toroidal field 
at a latitude of 60° versus the diffusivity at the convection layer for different values of the tachocline width has revealed that these fields are 
mainly eliminated for tachoclines with width t/[>0.08i? o (for rj c -2 x 10 9 - 1 x 10'° cm 2 s -1 and r c =0.7l5R Q ); or f° r rfi^0.02^ o and almost any 
value of j] c in the appropriate solar range. For intermediate values of rfi^0.04i? Q - 0-06i? Q , strong toroidal fields should survive at high latitudes 
in the butterfly diagram and those values are therefore not suitable. We have built butterfly diagrams for both a thin and a thick tachocline that 
best match the observations. We have also found that a prolate tachocline is able to reproduce solar-like butterfly diagrams depending on the 
choice of appropriate diffusivity profiles and tachocline width range. 
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1. Introduction 

Over the last few years, with the increase of the computational 
power and the improvement in the observational techniques, 
there has been a substantial advance in the knowledge of our 
magnetic star. However there is still a considerable number 
of open questions regarding the magneto-hydrodynamic pro- 
cesses that govern the large scale solar magnetic phenomena, 
i.e., the 11-years sunspot cycle, the polarity inversions of the 
magnetic field, the phase lag between the toroidal and poloidal 
inversions, and the latitudinal distribution of the toroidal mag- 
netic field. In the currently accepted scenario, the large scale 
solar magnetic cycle (LSSMC) is governed by a dynamo action 
that is the responsible for the transformation of a positively ori- 
ented poloidal magnetic field in a negatively oriented toroidal 
field, and the subsequent transformation of this one in another 
poloidal field but with the opposed polarity, and so on, until 
completing the cycle. The first stage of the process is the well 
known Q effect: the poloidal field lines are dragged ("frozen 
in"), and amplified by the differential rotation of the fluid. From 
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early results of the helioseismology, there has been a common 
agreement that this process takes place in the tachocline, a thin 
layer (< 5% of the solar radii) located at the base of the solar 
convection zone (~ QJR Q ) where the transition from a uni- 
formly rotating regime of the radiative core to a differentially 
rotating regime in the convective envelope occurs. The strong 
radial shear (dCl/dr) that exists in this region suggests that O 
acts mainly in the tachocline instead of in the entire convection 
zone. 

The second stage of the process, conventionally called a effect 
from the early days of the Parker turbulent mean field dynamo, 
consists in to convert this belt of toroidal field in a new poloidal 
field with opposite polarity. This subject has been the center of 
intense discussions and debates and an excellent revision of 
both, the a effect mechanisms and the s everal dynamo mod- 
els can be found in Charbonneau (2005) and references there 
in. In this work we assume an a effect as the result of the de- 
cay of active bipol ar magnetic reg ions (BMR), as it was orig- 
inally proposed by Bab cocld(ll96ll) and followed bv lLeightonl 
(1969). The fundamental ingredient in this mechanism is the 
emergence across the convective bulk of magnetic flux tubes 
due to the Parker-Rayleigh-Taylor instability in order to form 
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sunspots or active BMRs at the surface, while the BMRs de- 
cay they migrate to the poles carrying with them the necessary 
magnetic flux to neutralize the remnant poloidal field and gen- 
erate a new one. Magnetic reconnection between big loops of 

both hemispheres happens in this phase, as observed. 

Magnetic flux tube s imulations (D'Silva&Choudhuri 




Fan et al. 


1993; 


ICaligarietal1l995. 1998; 


Fan & Fisher 



5 x 10 4 and 10 5 G are able to erupt across the convection zone 
and emerge to the surface with the observed inclinations in the 
sunspots (joy's law). A sub-adiabatic layer is necessary in or- 
der to store magnetic field up to this strength and once again 
the tachocline is the best place to allow this storage and ampli- 
fication. The decay of the BMRs at the surface requires both, 
super-granular diffusion and meridional transport, for this rea- 
son the numerical models of the solar dynamo that work in 
the kinematic regime, that is, ignoring dynamical back-reaction 
of the magnetic field on the flux, include four fundamental in- 
gredients: differential rotation, meridional circulation, diffusion 
terms and a source of poloidal field (the a term). 

Recent kinematic models in the above scenario have 
been able to reproduce the majority of the LSSMC features, 
but have failed at producing a correct latitudinal distribu- 
tion of the toroidal field, showing intense magnetic fluxes at 
higher latitud es and, conseq uently, sunspots close to the poles. 
iNandv & Choudhuril j2002l) found a possible solution to this 
problem by allowing the meridional flow to penetrate below 
the tachocline. Under this hypothesis, the magnetic flux will 
be stored in a highly sub-adiabatic region and will emerge 
only at the desired latitudes. This assumption has given rise 
to a new controversy: how much of the meridional flow must 
penetrate below the tac holine? On one side, in a recent work 
IChatteriee et al . (2004), working under this hypothesis, have 
been able to reproduce the characteristics of the solar cycle in 
the two hemispheres, and also the observed parity rule (Hale's 
law). On the other side, numerical simulati ons of a m eridional 
flow penetrating in a sub-adiabatic medium (Gilman & Miesch 
2004) have shown that the dynamical effects of the fluid alone, 
without including magnetic field, do not allow a penetration 
below 5% of the tachocline and, in the case that the flow 
can penetrate in some way the radiative zone, the penetra- 
tion would be reduce d to a bout f ew kilometers only. In a 
recent report, iRiidiger et alJ J2005I) confirm this result argu- 
ing further that a meridional circulation confined within the 
convection zone alone would be able to produce an aQ. dy- 
namo. Another problem that arises when the penetration of the 
flow in the radiative zone is allowed is the excessive burn- 
ing of light elements in such a hot zone. Numerical mod- 
els of mixing, based on helioseismic measurements of the 
sound speed and density profiles, indi cate a maximum extent 
of mixing of 5% ins ide the tachocline (Brun et al. 2002). Also, 
IGuerrero & Munozl J2004I) develo ped a hybrid model using the 
profile of meridional velocities of Nandv & Choudhuri (2002) 
and the a-term of Dikpati & Charbonneaul & 19991) and found 
that a deep meridional flow solves only partially the problem 
of the distrib ution of the toro idal fields at the polar regions. 
Recently, Di kpati et alJ J2004I) have found an appropriate com- 



servations, however, until now there is no clear physical pro- 
cess that is able to explain why the sunspots appear only at 
lower latitudes. 

In this work, we us e a modified version of t he numer- 
ical code developed by IGuerrero & Muno'zl (|2()04) introduc- 
ing more det a iled diffusivity and a profiles, as suggested by 
Dik pati et alJ ( 120041) . aiming to explore how the shape and 
thickness of the tachocline may influence the latitudinal distri- 
bution of the toroidal field. Although during the entire dynamo 
process there are several mechanisms that are still not under- 
stood, such as, the exact values and the variations of the mag- 
netic diffusivity in the convective and radiative envelopes, and 
the return flow properties, the actual radial and latitudinal oper- 
ation of the a-effect, or the quantity of poloidal field that can be 
produced by the decay of the BMRs, for our propose here, we 
will consider profiles that are constrained by helio seismology 
observations and also by new results of simulations. 

In section §2 we present the basic mathematical formal- 
ism of the solar dynamo mechanism and some details about the 
construction of the mo del (a complete description of the nu- 
merical tool is given in IGuerrero & Munoz (2004)); in §3, we 
present a detailed discussion of the assumed profiles and the set 
of parameters that allow for a best fit to the observations. The 
results of our analysis are described in §4; and finally in §5 we 
sketch our conclusions. 



2. The model 

The equation that describes the spatial and temporal evolution 
of the magnetic field in an ionized medium is the magnetic 
induction equation, which under the assumption of azimuthal 
symmetry can be divided in its poloidal and toroidal compo- 
nents, respectively: 



OA 1 



1 



— + -(U ■ V)(M) = // 1 V 2 - - | A + S , (r, 8, B) , 
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(Bp ■ V)Q 
-V?7X(V xB) + rj\V L - — \B , 



(1) 

(2) 



bination of parameters that are able to better reproduce the ob- poloidal field. 



where R p -B r + Bg-V xA is the poloidal field and B the toroidal 
field, Q is the angular velocity, u-u r +ug denotes the meridional 
components of the velocity field, s-r sin 9, r\ is the magnetic 
diffusivity, and S i (r, 9, B) is a source that describes the a-effect 
in the Babcock-Leighton dynamo. 

The above equations are solved in a 2-dimensional mesh of 
128 x 128 grid points by using the ADI implicit method within 
the physical domain: 0.55fl G < r < R Q , and 0° (pole) < 8 < 90° 
(equator). The boundary conditions are as follows: B and A are 
set to zero in the pole and in the bottom radial boundary. At the 
top boundary we set a vacuum condition, so that, B(R Q , 8) = 
and A(R Q , 8) is smoothly matched with a vacuum photosphere- 
corona, (V 2 - p s ) n 2 fl )^-0- At the equator, we set the toroidal 
component, B-Q, and a continuous boundary, §^=0, for the 
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Fig. 1. Isorotation lines of the solar interior inferred from he- 
lioseismological observations. The dashed lines show the ap- 
proximate thickness of the tachocline. 

3. The dynamo ingredients 

As discussed above, a kinematic solar dynamo model needs 
four fundamental ingredients: differential rotation, meridional 
circulation, magnetic diffusivity and a poloidal source term. 
The characteristics of these ingredients are in part constrained 
by helioseismology observations but there is some degree of 
freedom in the choice of the parameters that may reproduce 
better the LSSMC features. We assume the following profiles 
and parameters in order to generate a fiducial model able to 
reproduce the observed butterfly diagram. 

3.1. Differential rotation 

One of the goals of the helioseismology was the determi- 
nation of an accurate profile for the angular velocity both 
at the surface and in deeper layers. The results revealed a 
radiative core that is rotating with uniform velocity. This 
changes to a convective bulk that is rotating differentially 
with a retrograde velocity with respect to the radiative inte- 
rior at higher latitudes and a pro-grade velocity at lower lati- 
tudes. The interface between these two regimes is a thin layer 
called tachocline. An analytical profile of this di fferential rota- 
tion was introduced by MacGregor & Charbonneau ( 1997) and 
used for the first time in a B abcock-Leighton type dynamo by 
Dik pati & Charbpn neau|(|l 999|). We employ her e the same ex- 
pression used by Dikpati & Charbonneau ( 1999): 



Q(r,0) = Q c + - 



1 +erf 2 



r-R c 



(O s (0)-Q c ) , 



(3) 



where Q c /27r=432.8 nHz is the uniform angular velocity of the 
radiative core, Q s (6)=Q. eij +a2 cos2 ® +a 4 c °s 4 & is the latitudinal 
differential rotation at the surface with Q. eq l2n-A6Q.l nHz be- 
ing the angular velocity at the equator, a-il1n--(yl.9 nHz and 
a4/27r=-67.13 nHz, erf(x) is an error function that confines the 
radial shear to a tachocline located in R c =0.1R o of thickness 
d\ . Figure 1 depicts the differential rotation profile in the solar 
interior for the core and convective layers. 



Until now there has been n o consensus about t he location 
and thickness of th e tachocline dCor bard et alj200l|). Some ob - 
servational works ( Ant ia et alll998t lcharbonne au et alJl 999). 
though of limited resolution, have found indications that this 
width may vary with latitude and suggested that the tachocline 
could have a prolate shape. We will show in section 4 how the 
change from a spherical to a prolate shape of the tachocline and 
its thickness influence the outputs of the model. 

3.2. Meridional circulation 

In the present stage of the observational techniques, there is 
a set of independent measurements of the meridional circu- 
lation that confirms a surface poleward flu x of about 20 m 
s -1 , which pers i sts un ti l a depth of 0.85 . R™ JHathawaviri 996 ; 
Hathaway et al. 1996; Latushko 1996; Snodarass & Dailev 
ll996HGiles et al.ll997llKomm et al.ll993HBraun & FanfT9 98). 
However, even helioseismology has been unsatisfactory to 
measure a deeper flow. Mass conservation predicts an equa- 
torward return flow and a unique convection cell per merid- 
io nal quadrant. Here, we use the same analytical prescriptio n 
of lNandv & ChoudhurilJ2002l) and lGuerrero & MufiozN2004 
which was previously introduced by iDikpati & C houdhuri 
dl994 . and lchoudhuri et al]dl995l) : 



p(r)u = V x [ifr(r, 0)e ] , 

where ^ is a stream function given by: 

7T(r-R p ) 



ifrr sin 8 — (r - R p )i//o sin 



(4) 



(5) 



_(Rb-R p ) 

x (l_ e -A« e )(l_^-/2)) 

v Jir-ro)/!:] 1 



and p(r) is a density profile for an adiabatic sphere with a spe- 
cific heat ratio y-5/3 (polytropic index m=1.5), thus: 



PW = C 



(t-H" 



(6) 



The value of i/'o and C are chosen in such a way that 
the amplitude of the meridional velocity, ug at middle lati- 
tudes is t/o=i/'o/C=25 m s _1 . The assumed value for the other 
five parameters are: j6i=6.06 x 10 9 cirT 1 , /?2=4.6 x 10 9 cm _1 , 
e=2.0000001, r =(R Q - R min )/4.0, and T=3.47 x 10 1() cm. R mi „ 
is the radius of the bottom of the convective boundary,and R p 
is the radius of the penetration depth of the flow. 

In a previous work, it has been shown that a deep merid- 
ional flow solves only partially the prob lem of the distribu - 
tion of the toroidal magnetic field ( Guerrero & Munozl f2004). 
Besides, as stressed in §1, this assumption seems to need 
physical support and does not seem to be confirmed by nu- 
merical simulations either, so that, in the present analysis we 
will not adopt a de e p mer idion al distribution (s e e how ever 
iNandv & Choudhurll iho02l) . and IChatteriee etall d2004 . for 
an alternative interpretation). 

Th e dynamical shallow water model of Gilma n~& Mieschl 
d2004 showed that it is impossible for the meridional flow 
to go below a few percent of the top for the solar tachocline. 
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Fig. 2. The adopted meridional circulation lines profile in this 
work. The dashed lines represent the tachocline with a constant 
thickness di=0.05R Q in this case. 

Recent calculations of Riidi ger et alJ J2005) support this result 
arguing that a weak penetration is able to produce a flux dom- 
inated solar dynamo inside the solar convection zone, without 
the necessity of participation of the tachocline, however they 
did not consider the overshoot effect of plumes and jets that 
actually can penetrate inside the tachocline and even down in 
a thin fraction of a more sub-adiabatic medium jRoeers et alJ 
120061) . In our calculations, we will assume a weak penetration 
inside the tachocline and, as the width of the tachocline is al- 
lowed to vary in this work, the percentage of penetration will be 
able to vary, as well, but we will consider a constant value for 
the penetration radius R p -0.69R Q . Figure[2]shows the assumed 
meridional circulation profile. 



compatible with the values that have been obtained from he- 
lioseismic measurements, the fir st with a ve ry high precision 
and the sec ond as an upper limit dBasull997h . and r cl =O.95fl , 
d 3 =0.01R Q (Dikp ati etal1l2006l) (see Figure 0. Our model is 
specially sensitive to the the diffusivity variation between the 
radiative core and the convective envelope. With the choice 
of the parameters above, we are considering a weak turbulent 
regime for the lower region of the tachocline and a turbulent 
regime for the upper part of it, with a sharp variation between 
them, in order to allow both storage of magnetic field and over- 
shooting of material below the tachocline (see, however, §4.3 
for a more detailed analysis of the dependence of the model 
with the parameters r\ c and r c ). 
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Fig. 3. Adopted depth dependence of the magnetic diffusivity 
profile. Y axis is in log-scale and is normalized to the super- 
granular surface value (10 12 cm 2 s~'). 



3.3. Magnetic diffusivity 



The radial dependence of the magnetic diffusivity is proba- 
bly the most undetermined of the profiles of the solar inte- 
rior. Diffusion must exist in the convective envelope due to 
the intense turbulence present, but besides the values of super- 
granular diffusion observed at the surface (j]s ~ 10 12 - 10 14 cm 2 
s _1 ), the depth dependence is uncertain and, in general, only ap- 
proximate values are assum ed for this quantity i n the dynamo 
models. We will here follow D ikpati et al] f2004) and assume a 
weak turbulent diffusivity regime for the radiative interior, with 
r\ r -2.2 x 10 s cm 2 s _1 , a turbulent regime for the convective en- 
velope, with ?7 C =5 x 10 9 cm 2 s above the tachocline, and a 
third region of super-granular diffusivity, from r c \-Q.95R Q on, 
where rjs = 10 12 cm 2 s _1 . 



ij(r) = r] r + 



1c 



l+erf 



'Is 



1 + erf 



r-r c \ 
c/3 



■(7) 



Different combinations of the parameters d%, di,, r r , and r r \ 
in the e quation above can be considered (see, e.g. jDikpati et alJ 
The boundary between a weak and a strong turbu- 
lent regime is the location where an abrupt change in the tem- 
perature gradient occurs in helioseismic calibrations. We have 
presently considered r c =0.715/? G and fl?2=0.017? G which are 



3.4. The source of poloidal field (the a term) 

From the Parker mean field dynamo theory days until now, the 
a source term of the poloidal field remains virtually undeter- 
mined. As a coupling factor between eqs. ([0 and (0, the main 
function of this term is to generate poloidal field from pre- 
existent emergent magnetic flux tubes in the toroidal direction, 
and its form is oriented to resemble the Babcock-Leighton con- 
cept of formation of dipolar like field as the result of the decay 
of BMRs and the reconnection of loops from both hemispheres. 

The magnetic flux tubes, at the base of the convection zone 
are amplified by the differential rotation until they reach a value 
such that the magnetic pressure begins to dominate over the 
gravitational pressure. Numerical simulations suggest a value 
of 10 5 G. Then, they are lifted on the entire convection zone 
and twisted by the Coriolis force to form BMRs at the observed 
latitudes. 

In order to allow this process to occur, we let the poloidal 
field at the surface (r-R Q ) to be proportional to the toroidal 
field at the base of the solar convection zone (r-R c + 
d\ 12), where d\ is the t hickne ss of the tachocline (see, e.g. 
Dik pati & Charbonneaul Jl999l) V Since the toroidal field lines 
are mainly concentrated in a belt around the solar equator, 
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following Dikpati et alJ J2004I) . we consider the action of the 
Coriolis force also concentrated at lower latitudes, thus: 



Table 1. Values of the parameters used in the model of Figure 
5. 



S r (r,e,B) = a - 



1 +erf 



x sin 9 cos ( 



r — ri 
1 



1 - erf 



1 +e yi(n/4-ff) 



r~r 3 
B(R c ,0) 



(8) 



Bo 



2\ 



where r 2 =0.95fl G , r 3 =fi , d 2 =d 3 =0.01R Q , yi=30. The ampli- 
tude of the poloidal source is determined by ao for which we 
assume a fixed value of 130 cm s -1 . Figure|4]depicts the radial 
and latitudinal dependence of this profile. 

The latitudinal dependence of S i is of great importance in 
the reproduction of the features of the solar cycle. We tested al- 
ternative possibilities reported in liter ature, i.e., profiles which 
are proportio nal to sin 9 c os 9 (Dikpati & Charbonneau 1999), 
sin 2 9 cos 9 llKiiker et al.l 120011) . and cos 6 ( Chatt eriee et alJ 
2004) factors, but we found that the profile prescribed by eq. 
<03i is the one that better reproduces an observed butterfly dia- 
gram. 

The last term on the RHS of eq. ((HJ is a quenching term lim- 
iting the growth of the poloidal field. Actually, this quenching 
should be given by the back reaction of the magnetic field on 
the velocity field, but the physical way by which the poloidal 
field stops growing is unknown. This form, which was also as- 
sumed previously by many authors, simply reproduces the fact 
that a toroidal field with a value larger than Z?o=10 5 G would 
generate suns pot pairs with tilts which is in disagreement with 
the joy's law JP'Silva & Choudhurilll993l) . It is important to 
note that this term is the only source of non-linearity of the 
system and works non-locally. 

Other mode l s, like , e.g. iNandv & Choudhuril d2002l) and 
IChatteriee et al.l (|2004) employ a numerical procedure to in- 
clude the quenching of the poloidal field. In these models, 
whenever the toroidal magnetic exceeds Bo at the base of the 
solar convection zone, a fraction / of it is artificially made 
to erupt to the surface at time intervals r = 8.8xl0 5 s 
jNandv & Choudhurill200lh . In essence, this buoyancy mech- 
anism does the same as our quenching term above, mainly with 
a supergranular diffusion value at the surface, with a time delay 
in order to ensure the right phase relation between the toroidal 
and the poloidal fields, but this phase lag can be fitted by tuning 
the velocity of the meridional flow at the base of the convection 
zone. 

We notice that in a recent work, Dikp ati etal J f2004) have 
introduced an a effect located in the tachocline. It could occur 
as a consequence of hydrodynamic instabilities at the top of the 
tacholine and become important in two different ways. First, as 
this term does not depend on the toroidal field magnitude, it 
could be a source of field in times when activity goes down 
(i.e. during the Mounder minimum), second, it could provide 
the correct parity relation (Hale's law) when the integration of 
the equations is made in the two hemispheres. Both of these 
issues are out of the scope of this work and, for this reason 
have not been considered. 



PARAMETER 


VALUE 




2n x 460.7 nHz 


Rc 


O.7R 


di 


0.05fl o 


U 


25 m s" 1 


R P 


0.69fl o 




2.2 x 10 8 cm 2 s _J 


>h 


5.0 x 10" cm 2 s" 1 


n s 


1.0 x 10 12 cm 2 s" 1 


a 


130 cm s" 1 



4. Results 

Employing the profiles described above, we have first built the 
butterfly diagram of the model depicted in Figure [5] and Table 
1, with a constant width tachocline with d\ = 0.05/? Q . It re- 
produces some of the main features of the large scale solar 
dynamo model. That is, the periodicity of the magnetic cycle, 
the observed magnitude of the toroidal fields near the equator, 
as well as, the weak radial fields near the pole, and the phase 
lag between them. However, we notice that large toroidal fields 
persist above 45 degrees suggesting that sunspots should also 
appear at those latitudes, which is not observed. In the next 
section we explore alternative possibilities to this solution. 

4.1. An ellipsoidal tachocline 

As stressed in §1, so far, the exact location and width of 
the tachocline are still unknown. A good revision of the ob- 
servational results and theoretical mo dels of the tachocline 
can be found in ICorbard et alJ d200l|). Some evidenc e for 
a prolate shape has been found bv lAntia et ail dl998b and 
ICharbonneau et all Jl999l) . but this could be due to observa- 
tional uncertainties or to a particular sensitivit y to the inver- 
sion techniques employed to an alyze these data ijCorbar d et alJ 
2001). Numerical simulations dPikpati & Gilmaal200ll) have 
found that a strong magnetic stress (B > 10 5 G) can pile the 
matter into the pole increasing the density and thus the width 
of the tachocline there. This could explain, in principle, the pro- 
late shape, however, since a strong field is present only during 
the maximum of activity, this could suggest that the tachocline 
would have a varying shape through the cycle, becoming nearly 
prolate during the maximum activity. 

In order to evaluate the effects of a potential prolate 
tachocline in a flux dominated dynamo model, we have intro- 
duced a latitudinal dependence on the width of the tachocline 
in equation 0, making it to vary from fl?i(pole)=0.077? o to 
iii(equator)=O.O2/?0 (see left panel of Figure|6] top). This mod- 
ification in the shape of the tachocline has the effect of de- 
creasing the radial shear dCl/dr at the high latitudes. If the ap- 
pearance of strong toroidal fields at high latitudes were mainly 
sensitive to the radial shear, then one should expect that a re- 
duction in d£l/dr at those latitudes would reduce the ampli- 
fication of the toroidal field. However the butterfly diagram 
depicted in Figure |6] (top), which was calculated for a pro- 
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Fig. 4. Assumed radial and latitudinal profiles for the a effect mechanism. Radially, the a term is concentred close to the surface. 
Its latitudinal distribution corresponds to the belt where the sunspots appear. 
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Fig. 5. Time-latitude butterfly diagram for the model of Table 1 with a tachocline of constant width d\ =0.05/? Q . The continuous 
(dashed) lines represent the positive (negative) strength of the toroidal field at the base of the solar convection zone (i.e. the top 
of the tachocline r-R c + di/2). The lines are log spaced and cover the interval between 5 x 10 4 - 10 s G. The background gray 
scale represents the positive (clear) and negative (dark) radial field at the surface. 



late tachocline, shows no significant changes with respect to 
the diagram of Figure|5] for which a constant width tachocline 
was assumed. Besides, when an oblate configuration is consid- 
ered instead (left panel of Figure |6] bottom), we find that the 
toroidal field is amplified to its maximum value only below the 
60° latitude (right panel of Figure [6] bottom). In other words, 
for an oblate tachocline, for which the radial shear is improved 
towards the higher latitudes, we find an inhibition in the gen- 
eration of toroidal fields at those latitudes, contrary to what is 
commonly expected. Also, for both the oblate and the prolate 
tachoclines, the generation of toroidal fields at the low latitudes 
is practically the same, suggesting that <9Q /dr is not influencing 
the behavior of the toroidal field. 

4.2. A thinner or a thicker tachocline? 

Two interesting remarks can be pointed out from the results 
of Figure [6] On one hand, it does not seem to be posssible to 
reproduce the butterfly diagram with a prolate tachocline (i.e., 
with a thicker tachocline at higher latitudes), at least not for 
the assumed conditions. On the other hand, we find that the 
radial shear dtl/dr does not seem to be contributing for the 
amplification of the toroidal magnetic field. 



In order to investigate the origin of this apparent paradox, 
we have considered again a spherical tachocline (with constant 
width), and then computed the components of the shear term 
on the RHS of equation @. (B p • V)Q.=B r d£l/dr + B g /r dQ./d6, 
at a high latitude from the equator (0=60°) for r=R c (i.e., in 
the center of the tachocline where the radial shear is maximum) 
at a time when the radial magnetic field (B r ) inverts its polar- 
ity. At this time of the cycle, the action of the shear term upon 
the poloidal field generates new branches of toroidal field and 
establishes the morphology of this growing branch for the next 
phase of the cycle. We have also computed the toroidal mag- 
netic field (B), at the same latitude when it reaches its maxi- 
mum at the top of the tachocline. These quantities were com- 
puted as a function of the width of the tachocline (d\) taken in 
the range of possible values that are inferred from the observa- 
tions (0.01/?o <d x < 0.1j^: ICorbard et al.lfeOOll) '). the results 
are depicted in Figure^ 

We note that the radial component, B r {^) (Figure^), ac- 
tually tends to decrease with the increase of the width, as ex- 
pected. But, more significant is the fact that its value is about 

1 This is the latitude above which strong toroidal fields > 5 x 10 4 
G should not appear, but have developed in the butterfly diagram of 
Figure[6](top). 
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Fig. 6. Angular velocity profiles and time-latitude butterfly diagrams for a prolate tachocline (top panel); and an oblate tachocline, 
(bottom panel). See the text for details. The contours specifications in the right panels are the same as in Figure[5] The tachocline 
in the left panels is represented by dashed lines. 



two orders of magnitude smaller than that of the latitudinal 
component ^r(^) (Figure^). Therefore, when a new toroidal 
field begins to be generated its growth is dominated by the lati- 
tudinal shear in eq. Tests run for different latitudes revealed 
the same effect. This is confirmed by Figure 7c, which shows 
that the magnitude of the generated toroidal field at high lati- 
tudes varies with the width of the tachocline in a similar way to 
the latitudinal shear component (Figure lb), attaining a maxi- 
mum value at a width, d \ that depends on the assumed magnetic 
diffusivity. This latter result will be explored in more detail in 
the next paragraphs. 

The results above raise two new questions. Since the ra- 
dial shear does not seem to be important for the amplification 
of the toroidal magnetic field, nor even at high latitudes, is 
the tachocline really participating in the dynamo process? And 
if yes, what is the real thickness of this layer? In the present 
model, the tachocline is not only the interface where the radial 
shear dtl/dr is maximized, but also the place where the merid- 
ional flow penetrates and the magnetic flux tubes are stored and 
amplified before erupting by buoyancy effects to the upper lay- 
ers, so that our answer to the first question is yes. In order to an- 
swer to the second question we can examine the amplitudes of 
the generated toroidal field at the latitude 60°, at the top of the 
tachocline in Figure 0:. It indicates that values of d\ <O.O2/? 
and di>0.08R Q 2 are appropriate to prevent the formation of 
strong toroidal fields at high latitudes, as required by the ob- 
servations. But which value of d\ is the best one? As we have 
mentioned above, the toroidal fields that develop at the high 
latitudes are also dependent of the adopted diffusivity profile. 
For a definitive answer, we have to wait for better helioseismic 
observations in the future, nonetheless, a more detailed scan- 



2 As it will be shown below, this result holds only for values of 
r] c <2x 10 10 cm 2 s _1 . 



ning of the diffusivity parameters, as shown in the next section, 
can shed some light on this question. 

4.3. Parameter dependence 

The present model employed a large number of parameters and 
we have performed some tests in order to explore how sensi- 
tive the results above are to them. The meridional circulation 
term in equations ©-(|6|l has two free parameters, the ampli- 
tude of the meridional flow at the surface, at middle latitudes, 
Uo, and the depth of the penetration of the meridional flow, 
R p . While the first does not affect the strength of the gener- 
ated magnetic fields (though it is very important to establish 
the period of the cycle), the second is able to change the lat- 
itude of formation o f stron g magnetic fields, as stressed by 
iNandv & Choudhuril d2002l) . We should note, however, that in 
this analysis we are assuming a weak penetration regime, and 
small variations around the value depicted in Table ^ do not 
change significantly the results above. 

The a term, whose semi-empirical profile was chosen in 
order to better reproduce the observations (Figure^, has only 
one free parameter, the amplitude, ao (in equation[8}. It deter- 
mines the exact amount of poloidal field to be regenerated in 
order to support the cycle. We find that our model is generally 
very insensitive to variations to it. 

For the diffusive terms of equation |7}, it is possible to es- 
tablish some constraints in the different regimes. For the most 
external layer, we have adopted a value for the diffusivity that 
has been obtained from the observations, so that it has been 
fixed in all the simulations. The value of the diffusivity at the 
radiative zone is not expected to affect the results either, but 
these can be sensitive to the assumed turbulent diffusivity at the 
convection zone (tj c ), as indicated by Figure^] We may choose 
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Fig. 7. Amplitude of the shear terms: a) fi,-(f?); b) ^ (^) at 
the time of inversion of the radial magnetic field, at a high lat- 
itude (60°), at the center of the tachocline and c) the toroidal 
field B at the time it reaches its maximum at the top of the 
tachocline as a function of the thickness of the tachocline. 
Dashed and continuous lines represent two different values of 
the turbulent diffusivity in the convection zone, rj c -5 x 10 9 cm 2 
s _1 , and 2 x 10 10 cm 2 s _1 respectively (the continuous line cor- 
responds to the value of Table [Q. 



an appropriate value for rj c based on the period of the cycle and 
the magnitude of the generated fields. If the value of tj c is too 
large, the system will enter in a diffusion dominated regime, 
thus reducing the period of the cycle in the butterfly diagrams. 
On the other hand, if rj c is too small, it will lead to too large 
values of the radial field at the poles. 

In Figure [3] we have plotted the maximum toroidal mag- 
netic field at the top of the tachocline, at a latitude of 60° (as in 
FigureEt), as a function of the diffusivity tj c , taken in the range 
of values which are appropriate for the solar cycle. Different 
line styles correspond to different values of the tachocline 
width. The dotted line at 5 x 10 4 G marks the limit between 
buoyant and no-buoyant magnetic flux tubes. If a curve is lo- 
cated above this limit, magnetic contours will appear above 60° 
in the butterfly diagram, and this is not desired. The top and 
bottom panels of Figure [3] have been plotted for two different 
values of the parameter r c of the turbulent diffusivity profile of 
the convection zone (eq.0. According to Figure [3J this is the 
transition radius between a weak and a strong turbulent regime. 
Though this value has been determined from helioseismic ob- 



servations with high precision (r c =0.715/? Q ), we have decided 
to vary it by 0.5% in order to check the sensitivity of the results 
to it. In the top panel, r c =0.715/? o is the same value adopted 
in the previous Figures. In the bottom panel, we have displaced 
this value to r c -Q.12R Q which means that a larger portion of the 
tachocline must lie in the less turbulent (sub-adiabatic) zone. 
Both, the top and the bottom panels indicate a similar behaviour 
for the different values of the tachocline width, but with a slight 
shift of the curves of the bottom panel towards larger diffusivity 
values, which are allowed only if the turbulent zone is smaller. 

The top panel shows that a tachocline with a width of about 
2% of the solar radius or less will produce solar like butterfly 
diagrams for almost the entire diffusivity range. Intermediate 
widths, d\ ^O.O4/? -O.O6/?0, are out of the allowed range of 
magnetic fields for any diffusivity, and larger widths between 
~O.O8/? to ~0. IR , will also produce butterfly diagrams which 
are in good agreement with the observations for r\ c from 2 x 10 9 
cm s~ 2 to 1 x 10 10 cm s~ 2 . The bottom panel also indicates 
that intermediate widths (afi^0.04R o -0.06/? o ) will produce in- 
appropriate butterfly diagrams, while a thin enough tachocline 
(di<0.02R Q ) will produce solar like butterfly diagrams for the 
entire range of appropriate diffusivities, and thicker tachoclines 
d\ -O.O87?0-O. 1R Q will produce solar like results only for diffu- 
sivities in the 6xl0 9 -1.6xl0 10 cm 2 s _1 interval. Computations 
made with r c =0.1 IR G revealed a similar behaviour with a slight 
shift of the curves towards lower values of the diffusivity with 
respect to the r c =0.715 panel (Figure|8] top). 

The butterfly diagrams plotted in Figure|5]for two spherical 
tachoclines, one with a thin width (t/^0.02^0, 77 e =5 x 10 9 cm 2 
s -1 , r c =0.72R Q ) and the other with a thick width (di-0.08R e , 
T] c -3 x 10 9 cm 2 s , r c =0. 715,/?0) are the ones that better repro- 
duce the observations. 

We notice that the results above are also applicable to pro- 
late and oblate tachoclines and naturally explain the results of 
Figure |6] In fact, since a prolate tachocline has a larger width 
at higher latitudes (top panel of Figure |6), then the toroidal 
field contours with strength between 5 x 10 4 G and 1 x 10 5 
G develop over the entire hemisphere because the latitudinal 
shear (which dominates over the radial shear in all latitudes) 
increases towards the poles. On the other hand, in the oblate 
case (bottom panel of Figure [6} where the width is smaller at 
higher latitudes, the toroidal fields are suppressed there due 
to the smaller latitudinal shear, thus resulting in a more con- 
centrated field distribution at the lower latitudes. However, the 
thickness of the tachocline at which the latitudinal shear has 
a maximum or minimum value depends on the assumed mag- 
netic diffusivity profile. For the one assumed in Figure|6] a min- 
imum value is obtained for a thinner tachocline (as in Figure|9] 
top) and therefore, the oblate configuration is the one that better 
reproduces the observations, but this scenario could change if 
a different diffusivity value had been adapted for the tachocline 
zone, than the one assumed in Figure [6] Note that the param- 
eters used to built the butterfly diagram of the prolate config- 
uration (d[=0.07R Q at the poles, tj c =5 x 10 9 cm 2 s _1 ) lie in 
the forbidden zone of Figure [5] For example, according to the 
top panel of Figure [3] if we had taken tj c between 2 x 10 9 cm 2 
s _1 and 1 x 10 10 cm 2 s _1 , then a prolate configuration with a 
tachocline width at the poles with d\ between (0.08 - Q.l)R Q 
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Fig. 8. Maximum of the toroidal magnetic field at the top of the 
tacholine as function of the diffusivity (in log-scale)at a latitude 
of 60°, The different line styles represent different widths of 
the tachocline d\ . The doted line represents the limit between 
buoyant and no-buoyant magnetic fields 5 x 10 4 G, as explained 
in the text. Only the values below this line will appear in the 
desired latitudes. The top and bottom panels correspond to two 
different values of r c in the diffusivity profile, r c -Q.115R G (as 
in all previous Figures) and r c -0.12R Q , respectively. 



would produce an appropriate butterfly diagram, because in this 
case even at the high latitudes, where the tachocline is thicker, 
the latitudinal shear would be small enough to suppress the 
toroidal fields there (like in Figure |9j bottom). 

5. Conclusions 

In this work we have explored the effects of variations in 
both the shape and the width of the solar tachocline in a 
flux-dominated kinematic solar dynamo model. First, em- 
ploying an improved version of the numerical approach of 



iGuerrero & Muno'zl d2Q04l) with a choice of more realistic dif- 
fusion and a-effect profiles and assuming a tachocline with 
constant width d\ =O.O5/? , we were able to reproduce success- 
fully some of the main features of the 1 1 -year large scale solar 
magnetic cycle, like the phase lag between the toroidal and the 
poloidal fields, the correct period and magnetic field magni- 
tudes (Figure[5)- However remains of toroidal field component 
persisted at the high latitudes of the butterfly diagram for these 
initial conditions. 

Then, considering a prolate tachocline (Figure[6] top), with 
a larger width in the polar region, which implies a smaller ra- 
dial shear d£l/dr at these latitudes, we have obtained a butterfly 
diagram with similar distribution of the toroidal fields to that 
obtained with the constant width tachocline of Figure [5] On 
the other hand, when considering an oblate tachocline (Figure 
|6j bottom), we obtained a toroidal magnetic field with a latitu- 
dinal distribution that is in better agreement with the observa- 
tions, with an absence of toroidal fields at latitudes higher than 
60°. 

In view of these surprising results, we computed the 
toroidal field and the radial and latitudinal components of the 
shear term (eq. |2j that is responsible for the amplification of 
the toroidal field, (B p • V)£l=B r dQ./dr + B g /rdQ./d0 as a func- 
tion of the width (di) of the tachocline at a high latitude, at a 
position where d£l/dr should be maximum. The results suggest 
that the latitudinal component of the shear term dominates over 
the radial term for producing toroidal field amplification. 

We have also found that these results are very sensitive to 
the adopted diffusivity profile, specially in the inner convec- 
tion zone (which is characterized by the diffusivity r\ c and the 
radius r c of transition between a weak and a strong turbulent re- 
gion). A diagram of the toroidal field at a latitude of 60° versus 
the diffusivity at the convection layer for different values of the 
tachocline width has revealed that these fields are mainly elimi- 
nated for tachoclines with width d\ >0.08/? G and r] c -2x 10 9 - 1 x 
10 10 cm 2 s _1 , for r c =O.715i? ; and t] c =6 x 10 9 -1.6 x 10 10 cm 2 
s _1 for r c =0.72R & ; or for d\<0.02R Q and practically any value 
of rj c in the appropriate solar range. For intermediate values of 
d\ ^0.04 - 0.06/?o, strong toroidal fields should survive at high 
latitudes in the butterfly diagram and those values are therefore 
not suitable. The best fits to the observed butterfly diagram are 
shown in Figure [5] for a thin and a thick tachocline. We also 
conclude that a prolate tachocline can produce solar like re- 
sults depending on the choice of the diffusivity profile and the 
adopted range of the tachocline width. 

Finally, we should note that the poloidal magnetic field 
magnitudes are correctly reproduced except for a branch of 
strong radial fields migrating to the equator (see white features 
in the background grey scale of the butterfly diagrams of the 
figures). The formation of this branch is related to our non- 
local implementation of the a-term. We will explore alternative 
forms for this term in future work. 
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Fig. 9. Angular velocity profiles and time-latitude butterfly diagrams for a thin (di=0.02R Q , i] c =5 x 10 9 cm 2 s -1 , r c =0.12R Q ) 
tachocline (top panel); and a thick (d\=0.Q%Ro, r/ c =3 x 10 9 cm 2 s , r c =O.715/? ) tachocline, (bottom panel). See the text for 
details. The contours specifications in the right panels are the same as in Figure[5] The tachocline in the left panels is represented 
by dashed lines. 



References 

Antia, H. M, Basu, S., & Chitre, S. M. 1998, MNRAS, 298, 
543 

Babcock, H. W. 1961, ApJ, 133, 572 
Basu, S. 1997, MNRAS, 288, 572 
Braun, D. C. & Fan, Y. 1998, ApJ, 508, L105 
Brun, A. S., Antia, H. M., Chitre, S. M., & Zahn, J.-P. 2002, 
A&A, 391,725 

Caligari, R, Moreno-Insertis, F., & Schussler, M. 1995, ApJ, 
441, 886 

Caligari, R, Moreno-Insertis, F, & Schussler, M. 1998, ApJ, 
502, 481 

Charbonneau, R 2005, Living Rev. Solar Phys., 2, 2 
Charbonneau, P., Christensen-Dalsgaard, J., Hening, R., et al. 

1999, ApJ, 527, 445 
Chatterjee, P., Nandy, D., & Choudhuri, A. R. 2004, A&A, 427, 

1019 

Choudhuri, A. R., Schussler, M., & Dikpati, M. 1995, A&A, 
303, L29+ 

Corbard, T., Jimenez-Reyes, S. J., Tomczyk, S., Dikpati, M., 
& Gilman, P. 2001, in ESA SP-464: SOHO 10/GONG 2000 
Workshop: Helio- and Asteroseismology at the Dawn of the 
Millennium, ed. A. Wilson & P. L. Palle, 265-272 
Dikpati, M. & Charbonneau, P. 1999, ApJ, 518, 508 
Dikpati, M. & Choudhuri, R. A. 1994, A&A, 291, 975 
Dikpati, M., de Toma, G, Gilman, P. A., Arge, C. N., & White, 

O. R. 2004, ApJ, 601, 1136 
Dikpati, M. & Gilman, P. A. 2001, ApJ, 552, 348 
Dikpati, M., Gilman, P. A., & MacGregor, K. B. 2006, ApJ, 
638, 564 

D'Silva, S. & Choudhuri, A. R. 1993, A&A, 272, 621 
Fan, Y. 2004, Living Rev. Solar Phys., 1, 1 



Fan, Y. & Fisher, G. H. 1996, Sol. Phys., 166, 17 

Fan, Y, Fisher, G. H., & Deluca, E. E. 1993, ApJ, 405, 390 

Giles, P. M., Duvall, Jr., T. L., Scherrer, P. H., & Bogart, R. S. 

1997, Nature, 390, 52 
Gilman, P. A. & Miesch, M. S. 2004, ApJ, 61 1, 568 
Guerrero, G. A. & Munoz, J. D. 2004, MNRAS, 350, 317 
Hathaway, D., Gilman, P., Harvey, J. W., et al. 1996, Science, 

272, 1306 
Hathaway, D. H. 1996, ApJ, 460, 1027 

Komm, R. W., Howard, R. F, & Harvey, J. W. 1993, Sol. Phys., 
147, 207 

Kiiker, M., Rudiger, G, & Schultz, M. 2001, A&A, 374, 301 
Latushko, S. 1996, Sol. Phys., 163, 241 
Leighton, R. B. 1969, ApJ, 156, 1 

MacGregor, K. B. & Charbonneau, P. 1997, ApJ, 486, 484 
Nandy, D. . & Choudhuri, A. R. 2001, ApJ, 551, 576 
Nandy, D. . & Choudhuri, A. R. 2002, Science, 296, 1671 
Rogers, T. M., Glatzmaier, G. A., & Jones, C. A. 2006, astro- 
ph/0601668vl 

Rudiger, G, Kitchatinov, L. L., & Arlt, R. 2005, A&A, 444, 
L53 

Snodgrass, H. B. & Dailey, S. B. 1996, Sol. Phys., 163, 21 



